IUHET 243 
LA-UR-93-206 
January 1993 



SUPERSYMMETRY AND A TIME-DEPENDENT LANDAU SYSTEM 



V. Alan Kostelecky 

Physics Department 
Indiana University 
Bloommgton, IN 47405, U.S.A. 

V. I. Man'ko 

P. N. Lebedev Physical Institute 
Russian Academy of Science 
Moscow 111924, Russian Federation. 

Michael Martin Nieto 

Theoretical Division 
Los Alamos National Laboratory 

University of California 
Los Alamos, NM 87545, U.S.A. 

D. Rodney Truax 

Department of Chemistry 
University of Calgary 
Calgary, Alberta, Canada T2N 1N4- 

A general technique is outlined for investigating supersymmetry prop- 
erties of a charged spin-| quantum particle in time- varying electromag- 
netic fields. The case of a time-varying uniform magnetic induction 
is examined and shown to provide a physical realization of a super- 
symmetric quantum-mechanical system. Group-theoretic methods are 
used to factorize the relevant Schrodinger equations and obtain eigen- 
solutions. The supercoherent states for this system are constructed. 



1 Introduction 

The quantum behavior of a nonrelativistic charged spin-| particle in the presence of 
a constant and uniform magnetic induction is of importance in many physical con- 
texts. The wave functions for this system are solutions of the Pauli equation, which 
has two components corresponding to the two possible orientations of the spin. Each 
component has an energy spectrum consisting of a tower of equally spaced levels, 
called Landau levels |] . The two sets of Landau levels are degenerate, except for the 
ground state. This system is known to provide a physical realization of supersym- 
metric quantum mechanics [0, [|. The supersymmetry generators act to reverse the 
particle spin, thereby mapping one tower of Landau levels into the other. 

The classical motion of a point charge in a constant and uniform magnetic in- 
duction is rotation about a circle in the plane perpendicular to the magnetic field. 
This motion is most closely reproduced in the quantum system by coherent states 
j|. [5|, |6|, for which the expectation values of the charge's coordinates follow the clas- 
sical cyclotron motion. Coherent states can also be introduced for the spin-^ Landau 
system. The presence of the supersymmetry makes possible an extension of these 
states, resulting in supercoherent states [0]. 

It is natural to ask whether the notions of supersymmetry and of supercoherent 
states can be introduced in the context of the motion of a charged spin-| particle 
in more general electromagnetic fields. As the construction presented in Ref. M 
relies on the factorization of the hamiltonian, it is not apparent a priori how to 
handle more complicated situations. The present paper addresses this issue. We 
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demonstrate that a group-theoretic analysis can provide the key to a supersymmetric 
factorization. Here, we focus primarily on the case of a uniform but time-dependent 
magnetic induction as an explicit example. However, the formulation of the problem 
and the methods used are applicable in a broader context. 

We also use the results to obtain supercoherent states. Our construction ex- 
tends the previously developed coherent states for a spinless charge in time-dependent 
magnetic (and electric) fields || ||, [10| . This earlier approach used time-dependent 
integrals of the motion satisfying an oscillator algebra and the standard displacement- 
operator method [11, 12, 1~3| . 

In Sec. 2, we establish our notation and perform a first separation of variables, 
using the invariance of the Schrodinger operator under translations along the direction 
of the magnetic induction. A rotated variable set is introduced that simplifies much 
of the subsequent analysis. The group-theoretic analysis of the resulting equations is 
presented in Sec. 3. We seek symmetries of the problem using the methods detailed 



in Ref. |14| and applied in Refs. (15, 16), \L% . The complexified symmetry algebra 
is constructed in Sec. 3.2. This generalizes the dynamical symmetry group for the 



constant-induction case, presented in Ref. [18 



We use these results to develop several factorization schemes, which are given 
in Sec. 4. The solution to the relevant Schrodinger equation is obtained in Sec. 
5, using group-theoretic techniques and some representation theory taken from Ref. 
| I9| . In Sec. 6, we extend these expressions to solutions of the Pauli equation, and in 
Sec. 7 the supersymmetry is explicitly identified. Finally, we construct the relevant 
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supercoherent states in Sec. 8, thereby completing the generalization of the problem 
treated in Ref. JF]. The coherent states for a single tower of levels, allowing for the 
time variation, are contained as a limit of these supercoherent states and are closely 



related to those constructed from integrals of motion in Ref. ||20|| . The Appendix 
demonstrates the reduction of all these results to the time-independent case. 
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2 The Time-Dependent Landau Problem 

Consider a nonrelativistic spin-| particle of mass M and charge e moving with mo- 
mentum p in a time- dependent electromagnetic field with four- vector potential (0, A) . 
The Pauli equation for this system is 

{M° ■ (P - eA ( r > 0)P + e <K r > *)}*3(r, t) = idt* 3 (T, t), (1) 

where ^3 is a two-component wavefunction in three space dimensions and the quantity 
a = (a x , <Ty,a z ) is a vector consisting of the three Pauli matrices. 

In this paper, we shall assume that the scalar potential 0(r, t) is zero and that 
the vector potential A(r, t) describes a uniform, time- dependent magnetic induction 
B. For convenience, we work in in vacuo where the magnetic induction B is related 
to the magnetic field H by B = /U H, and we choose the cylindrical gauge 

A x = -\By, A y = \Bx. (2) 

Writing the upper and lower components of ^3 as and \& 3 _, the Pauli equation 
reduces to the two equations 

{T 2 X + T 2 y + Tl + 2Me0 =p eB - 2iMd t )V 3± = 0, (3) 

where 

T x = p x + eBy/2, T y = p y - eBx/2, T z = p z . (4) 

Substituting for p x , p y , and p z the usual operator forms, these equations can be 
rewritten as 

S 3± ^ 3± = 0, (5) 
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where 

S 3± = d xx + d yy + d zz + iwL z - \w 2 {x 2 + y 2 ) T w + 2id r (6) 
are called the Schrodinger operators in three space dimensions and where we have set 

H = 1 , r — t/M , eB{r) = w(r) , L z = yd x - xd y . (7) 



The operators S$± commute with the ^-translation operator d z . This implies [13 
that we can separate Eq. (|5]) with respect to z. Set 

*3±(r,r) = V±{x,y,T)Z K± {z) . (8) 

The functions Z K± (z) satisfy the eigenvalue problems 

- id z Z K± = k±Z k± , (9) 

with the usual plane wave solutions. This procedure reduces Eq. (^) to 

S 2 ±^± = , (10) 

where the Schrodinger operators i>2± m t wo space dimensions are 

S 2 ± = d xx + d yy + iwL z + 2id T - 2h 2 (x 2 + y 2 ) - 2h ±, (11) 

with 2h ± = ^pw + k 2 ± and 2/z 2 = w 2 /A. These new equations depend only on the 
variables x, y, and r. The quantities k± are the constants of separation. 

It is convenient for our analysis to eliminate the iwL z term from (pp . Introduce 
the operator 

R=exp[r ) (r)L z ], (12) 



where rj (r) is to be chosen below to eliminate iw(r)L z from the expressions for S 2 ±. 
Define rotated solutions Q± by 

V ± = R- 1 Q ± . (13) 

The rotated Schrodinger operators S 2 ± are given by 

S 2 ± = RS 2 ±R~ 1 = exp[r)L z ]S 2 ±exp[—rjL z ] 

= d xx + d yy + iw(T)L z + 2id T -2tr]L z -2h 2 (x 2 + y 2 )-2h ±. (14) 

Setting 

V( T ) = I / w(fi)dfi (15) 

simplifies the expressions to 

S 2± = d xx + d yy + 2id T - 2h 2 {x 2 + y 2 ) - 2h 0± . (16) 

These operators are the Schrodinger operators for time-dependent isotropic harmonic 
oscillators. In the rotated frame, Eq. fliCf ) becomes 

5 2± 6 ± = 0. (17) 

Throughout the remainder of this paper, we use the usual roman letters to repre- 
sent operators in the original space of the problem (|i~0D , and script letters to represent 
operators associated with the rotated equations (|17D . 
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3 Lie Symmetries 

If w = eB is time- independent, the solution of Eq. (O) is possible by a direct treat- 
ment. In section IV of 0, the separation of variables was performed by setting p z = 
and ^±(x, y, r) = if)±(x, y)e~ lE±T ', and the resulting differential equation was factored 
into raising and lowering operators. However, when w varies with time it is no longer 
immediately apparent how to separate variables or how to identify appropriate raising 



and lowering operators into which the the differential equation (1Q|) can be factored. 
Instead, we proceed with a systematic approach that makes use of symmetries of the 
Schrodinger operators iS>2± in Eq. fllEl) . 

To simplify the expressions in this section and in Sees. 4 and 5, we write equations 
only for the case associated with the Schrodinger operator iS 2 + ; rather than both 5 2 ± 
at once. We also denote £2+ by S2 , and h 0+ by ha. Analogous expressions for the 
problem with <S>2_ can be found by replacing all occurrences of ho = h 0+ by /to-- 
This means that all the results obtained here and in the subsequent two sections have 
duplicate forms. We discuss the role and significance of this duality beginning in Sec. 
6. 

3.1 The Symmetry Operators 



The symmetries we seek for the Schrodinger operator S2 have the form |L4|, [15 
t = A(x, y, r)d T + B 1 (x, y, r)d x + B 2 (x, y, r)d y + C(x, y, r). 
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These operators generate space-time transformations. For these space-time transfor- 
mations to be symmetries of flTEp, they must satisfy the commutator [|TJ], [TJ| relation 



[S 2 ,C]=A{x,y,r)S 2 , (19) 

where A(x,y,r) is some function of the space-time variables. In the unrotated space 
of Eq. fllOD , this equation takes the form 

[S 2 ,L] = \(x,y,r)S 2 , (20) 

where we have defined 

L = R^jCR = A(x, y, r)d T + B\x, y, r)d x + B 2 (x, y, r)d y + C(x, y, r) (21) 

and 

X(x,y 1 T) = R- 1 A(x 1 y 1 r)R. (22) 

We next proceed to establish and solve a set of differential equations that determine 
the explicit form of Eq. (JTS|). 

Substituting ([16]) for S 2 and ( |T8| ) for C in Eq. (|T9"D, we obtain a system of partial 
differential equations for the coefficients A, B 1 , B 2 , and C: 

A x = A y = , 2B\ = 2B 2 y = A , Bj + = 0, 
B^ + BJ„ + 2iB* + 2C, = , #L + B&, + + 2C7 y = 0. (23) 
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We solve (p3|) in the usual manner to obtain 



A = A(t) , A = A, 



B 1 = \Ax + f3 1 2 y + e 1 {r) , B 2 = -f3\x + \Ay + e 2 {r), 



C = -\A(x 2 + y 2 ) - ie x x - ie 2 y + e(r) 



where the r-dependent coefficients .4., e 1 , e 2 , and e satisfy: 



[4 + 8/i 2 4 + Ah 2 A = 0, 



e 1 + 2/^e 1 = 0, 



e 2 + 2h 2 e 2 = 0, 



ie - \A + k A + h A = 0. 



(24) 



(25) 
(26) 
(27) 
(28) 



Since Eqs. ( PB"D and ( |2"TD have the same form, they are satisfied by particular 
solutions Xi( r ) an d X2{ T ) that also have the same form. The coefficient of k\ is zero 



in (Pq), so the wronskian W(xi,X2) must be constant [|15| . We choose to scale the 
two solutions so that 



W(xi, X2) = X1X2 ~ X1X2 = 1. 



(29) 



The general solutions then have the form 



e l (r) = /3 n Xi(r)+/3 12 X 2 (r), 



e 2 (r) = P 21 Xi (t)+(3 22 X2(t). 



(30) 



It is known [[DJ that if Xi an d X2 are solutions of (E^), then 



<^i( r ) = (Xi) 2 ; ^(r) = (x 2 ) 2 , fsir) = 2x1X2, 



(31) 



are particular solutions of fl25|). The general solution is the linear combination 



A(r) = P'iMr)} + P 2 {Mr)} + P 3 Ws(r)}, (32) 

where j3 2 , and (3 3 are real constants. 

At this stage, Eq. (EH) can be integrated to yield 

e(r) = /^{f^i + ^0^1} + /9 2 {i0 2 + ik oV2 ] + /5 3 {i0 3 + ^0^3} + /3 4 «- (33) 
The remaining coefficients in the Lie derivative are then found to be 

B 2 = -/3 2 H4+/3H|^} + ^{|^} + /5 3 {|^} + ^ 21 {Xi} + ^ 2 {X2} ) 
C = ^{-^(x 2 + y 2 ) + \^ 1 + ih oVl } + (3 2 {-{(p 2 (x 2 + y 2 ) + ±<p 2 + ih <p 2 } 
+ P 3 {-iV3(% 2 + y 2 ) + + iho<p 3 } + (3 u {-ix X i} + /9 12 {-^x 2 } 

+ P 21 {-iyxi] + P 22 {^yx2} + ^{i}- (34) 

Finally, the generators of the symmetry group of the Schrodinger operators S 2 can 
be obtained, by direct substitution. Three of them have the form 



Cj = ipjd T + \(pj{xd x + yd y ) - \(f>j{x 2 + y 2 ) + \<pj + ih Q ipj 



(35) 



for j = 1,2,3. These operators satisfy the commutation relations of an si (2, It) Lie 



algebra |15 | 



[£ 3 ,£ 1 ] = -2£ 1 , [£ 3 ,£ 2 } = 2£ 2 , [£ u £ 2 \ = £ 3 . 



(36) 



Another generator is L z , which spans an o(2) algebra and commutes with the £f 



[L Z ,£ J ] = 0. j = 1,2,3. 



(37) 
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The remaining five generators span a Heisenberg-Weyl algebra, 102, in two-dimensional 
space |T3|. These operators have the forms 



E = i, 

Ji = Xid x - ixxi, J2 = X2d x - ixx2, 

£1 = Xidy - iyxi, ^2 = X2d y - iyX2- (38) 
Their nonzero commutation relations are 

[Ji,Ji] = -£, [!C 1 ,!C 2 ] = -£. (39) 

The full symmetry algebra is the Schrodinger algebra in two space dimensions, 
(si (2, R) © o(2)) o w 2 . The remaining commutation relations are 

[£ h ji) = 0, [C2,Ji] = -J 2 , [C 3 ,Ji] = -Ji, 

[C 1 ,J 2 ] = Ji, [C 2 ,J 2 } = 0, [£ 3 , J%\ = J 2 , 
[£i,K,i) = 0, [C 2 , JCt] — -K 2 , [Cs,Kli] = — K,\, 

[C-i, ic 2 ] = iCi, [c 2 , jc 2 ] = 0, [£ 3 , ic 2 ] = ic 2 , 

[L z , J a ] = 1C a , [L z , K a ] = J a , a = 1,2. (40) 
3.2 Complexification of the Symmetry Algebra 

To work with Hermitian or Hermitian-conjugate operators, we need to complexify 
the symmetry algebra jlJJ. We begin by looking at the solutions of the differential 
equations (B^) and (^71). The real solutions to these equations were denoted Xi( T ) 
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and X2( r )- We can obtain complex solutions from these by defining 

£(r) = ^(Xi(r) + %(r)), £(t) = ^(Xi(r) - %(r)). (41) 



In this case, the general solutions to Eqs. and (]27|) can be written 

^(r) = /3 n {0 + /? 12 {a e 2 (r) = /3 21 {£} + (3 22 {£}. (42) 
The Wronskian of these solutions is 

W(t,£)=ti-it = -iW( X i,X2) = (43) 

Complex solutions to the differential equation (|25|) can be written in analogy to 
the real solutions (1317). We have 



Vi=t\ = ¥> 3 = 2&. (44) 

Then, in terms of these complex solutions, three of the generators are 

tj = ¥jd T + \<pj{xd x + yd y ) - \(pj{x 2 + y 2 ) + \<pj + ih ipj, (45) 

where j = 1, 2, 3. We can express these operators as linear combinations of the original 
operators (|3~5"D , which form the basis of the si (2, R) subalgebra, in the following way: 



A = |(A-£ 2 + iC 3 ) , C 2 = l(C 1 -C 2 -iC 3 ) , £ 3 = £ 1 + £ 2 . (46) 
It is more convenient to define the operators 

M 3 = iC 3 , M+ = £ 2 , M- = -Ci. (47) 
These operators satisfy the commutation relations 

[M+,M-] = -M 3 , [M 3 ,M+] = 2M+, [M 3 ,M-] = -2M— (48) 
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The operators fl4"T|) form a basis for an su(l, 1) algebra. 
The generator L z spans an o(2) algebra, as before. 
Finally, there are the five generators of a Heisenberg-Weyl algebra w%: 

J- = id x - ixi, J + = -ld x + ix£, 
K.- = id y - iyi, fC + = -£d x + iy£ } 

1=1. (49) 

These operators satisfy the nonzero commutation relations 

[J_,J + ]=I, [£_,£+]= I. (50) 

The operators J± and K± can also be expressed in terms of the operators J a and 
JC a , a = 1, 2, as follows: 

J- = -^{Ji + J+ = ^{-Ji + ij 2 ), 

£- = ^{£1 + 1X2), £+ = ^(-£i + i£ 2 ). (51) 
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The remaining commutation relations of the full Schrodinger algebra are 

[L z ,M 3 ] = [L z ,M + ] = [L z ,M_] = 0, 

[L z , J_] = JC-, [L z , J + ] = /C+, 
[L Z ,K,J\ = -J-, [L Z ,K, + ] = -J+, 
[M-,J-]=0, [M+,J-]=-iJ+, [M 3 ,J-] = -J-, 
[M-,J+] = -iJ-, [M+,J+] = 0, [M 3 ,J + ] = J + , 
[M-,1C-] = 0, [M+, K-] = -iJC + , [M 3 , K-\ = -XL, 
[M-,K+\ = -iK-, [M+,JC + } = 0, [M^K + \=K + . (52) 

It is also useful to have the expressions for the operators and equations in the 
original, unrotated frame. The operators of the Schrodinger algebra for (|T0|) can be 
obtained from fl2~l~|) and are given by 

j_ = R- X J_R = J_ cos{r]) - JC- wa.(ri), 
j+ = R~ X J+R = J + cos(?7) — /C + sin (77), 
= R~ 1 1C-R = KL- cos(r/) + J- sin(r/), 
k + = R~ l K, + R = /C+ cos(r]) + J+ sin(ry), (53) 
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and 

m 3 = R~ 1 M 3 R = M 3 + i<f 3 wL z 

= i{ip 3 {d T + \wL z ) + \ fa(xd x + ydy) - \(p 3 {x 2 + y 2 ) + ±y? 3 + ih ip 3 }, 

m + = R~ 1 M + R = M + + \(p 2 wL z 

= W2(d T + ~wL z ) + ~<p 2 {xd x + yd y ) - \(p 2 {x 2 + y 2 ) + \(p 2 + ih p 2 , 

m_ = R' l M^R = M- - \ifxwL z 

= -ipi(d T + \wL z ) - \(fi{xd x + ydy) + \(p\{x 2 + y 2 ) - \<p x - ih^. (54) 

The commutation relations (|47|) , (]50f) and ( |52D are preserved by the transformations 
(|5|) and (|§. 

A more convenient choice for a basis for W7 2 can be made. Define the operators 

a- = jsV- + ik.) = %{J. + iKJ) = + idy) - i(x + iy)t], 

c = j= 2 (j- - ik^) ^ e ^(J- - iK.) ^ ^[aO x - id y ) - i(x - iy)C). (55) 

and their conjugates a + and c+. These operators are related to the raising and 
lowering operators a, aJ, c, and introduced in Ref. J7J. Appendix A contains an 
explicit demonstration of the connection between the two sets of operators for the 
specific case of a charged particle moving in a constant magnetic induction. 

The nonzero commutation relations for the complexified Schrodinger algebra in 
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this new basis are: 

[m + ,m_] = — m 3 , [m 3 ,m + ] = 2ro + , [m 3 ,m_] = -2m_, (56) 

[a-,a+] = I, [c_,c+]=J, (57) 

[m 3 ,a±] = ±a±, [m 3 ,c±] = ±c±, (58) 

[m_, a + ] = — ic-, [m_, c+] = — ?a_, [m+, a_] = — ic+, [m+, c_] = — ia+, (59) 

[£*,a±] = =Fa ± , [£ z ,c±] = ±c±. (60) 
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4 Factorization of the Schrodinger Equation 

In this section, we continue to work with the Schrodinger operator S 2 = S 2 + as noted 
at the beginning of the previous section. Our next goal is to establish candidate 
factorizations for this operator. 

It suffices to work with the subalgebra Q = {m 3 , a±, c±, I}o{C z } of the Schrodinger 
algebra. If we denote the oscillator subalgebra generated by {m 3 , a±, c±, 1} as os(2), 
then Q = os{2) o o(2), where C z = iL z is the generator of o(2). The span of the 
oscillator subalgebra satisfies the nonzero commutation relations ([57]), (|58|), and (pOp . 

The first step is to calculate the operators a + a_ and c + c_. We find 



a + a_ = \[-\<fz{d xx + d yy ) + \<pz{xd x + yd y ) - i(yd x - xd y ) 

= \[-\<PsS2 + M3-C z -l\ = \[-\<pzS 2 + mz-C K -l\, (61) 



c+c_ = ~ [-\fi{d xx + dyy) + ~ip3(xd x + yd y ) + i(yd x - xd y ) 
+ ii(x 2 + y 2 ) + ^-l] 

= \[-\yzS 2 + Mz + C z -l] = |[-|^ 3 5 2 + m 3 + £*-!]. (62) 



and 



Rearranging Eqs. (^l|) and d62|) , we obtain 



- 2 ^S 2 



2a + a_ - M 3 + C z + 1 = 2c + c_ - .M 3 - C z + 1 



(63) 



-2V3S2 



2a + a 



m 3 + £ z + 1 = 2c + c_ — m 3 — £ z + 1 



(64) 
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where 

M 3 = i{(fsd r + \fa{xd x + yd y ) - {(p 3 {x 2 + y 2 ) + ±y? 3 + ih (p 3 }, (65) 

and m 3 is given by (|54|) . Equation (|63|) provides two different factorization schemes for 
the Schrodinger operator for the two-dimensional time-dependent harmonic oscillator, 
while Eq. fl64|) represents two factorization schemes for the Schrodinger operator for 
an electron in a time-dependent uniform magnetic induction. The equations are 
related to one another by the rotation R of (0). 

We remind the reader again that similar factorizations can be obtained for the 
Schrodinger operators S , 2 - an d «S 2 _. Note also that the operators a±, c± and L z are 
unaffected by the replacement of h 0+ by ho-. 
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5 Solution of the Time-dependent Schrodinger Equa- 
tion 

In this section, we present solutions of the two-dimensional Schrodinger equations 
obtained in Sec. 2. The results are found by group-theoretic techniques, and as such 
can be extended to more general situations than the one considered here. Once again, 
we remind the reader that, to minimize notational confusion, results are presented 
only for the operator S2 = S2+, the wavefunction \I/ = and rotated forms. The 
analogous expressions for £2- and can be found by replacing ho = h 0+ with h _. 

5.1 Further Separation 

Denote the solution space of the Schrodinger equation = by Qs 2 - The genera- 
tors of space-time transformations for the Schrodinger equation, say L of fl2~T|) , must 
satisfy Eq. (^). Then, for \I/ £ Q S2 , we see that 

[S 2 ,L]V = \S 2 ^ = 0. (66) 



This means the generators are constants of the motion . In particular, the gener- 
ators L z and 777.3 are two commuting constants of the motion. Therefore, we can find 
a set of common eigenfunctions, ^ u ,i, labeled by the eigenvalues u and I of m 3 and 
C Z1 respectively: 

m 3 ^ Uj , = u^ U)h C z ^ U)l = W U>1 . (67) 
Recall from (|T3| ) that \J/ Ui / = e~ vLz Q Ut i. Hence, we have 

M 3 e Uil = e^'mae-^'Q^ = uQ Ujl . (68) 
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Similarly, we obtain 

c z e u>l = ie Utl (69) 

for the second eigenvalue equation. 

In Eq. (|67f) , substitute the operator (|65|) , where 993 is given by (|4]). The resulting 
equation is a first-order partial differential equation for Q u ,i, which can be solved by 
the method of characteristics. The general solution has the form 

Q U:l = exp^^)> W(G, (i)TM, (70) 

where the separable coordinates are 



G = ^72, (2 = ^, » = T, (71) 

and where 



^=|03(Cl 2 + C 



I), 



T u ( / i) = ^3 1/2 (fr /2 exp(-^ ), 

Ldsh = K . (72) 



For details of the integration, see the Appendix in Ref. f|l|]]. Note that the function 



1Z cannot be written as a sum TZi(d) + 72. 2 (C2) + 7^3 (aO of arbitrary functions. The 



solutions Q U: i are called i?-separable solutions |TJ| of the equation S 2 Q u ,i = 0. 

Similarly, the functions that solve the first eigenvalue problem in (|67|) also 
solve the Schrodinger equation 5*2 ^ = 0. Since e~ vLz commutes with e %n , we obtain 
the solution to ([T0|) : 

*u,i = e^e-^^,/(Ci, (2)TM, (73) 
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where 1Z and T u are given by (|72|). Substituting for i in ( |I~0D and suppressing the 
u and / labels, we obtain 

fe: + fe 2 - (Ci + C2 )^ + 2# = 0, (74) 

which is the eigenvalue problem for a two-dimensional harmonic oscillator. 

A further separation of variables can now be made. We solve the eigenvalue 
equation (|69"D, where C z is given by 



C z = i{yd x - xd y ) = z(C 2 d Cl - d%), (75) 
and O u j by ([701). Since £ z commutes with 71, we obtain the eigenvalue problem 

£ 2 VV = tyu,h (76) 
which yields directly the first-order partial differential equation 

C2 dcr Ci wr-^- (77) 

Solving this equation by the method of characteristics, we get the solution 

ipu,l = T n ,/(p)e iW , (78) 
where the variables of separation p and 9 form a polar coordinate system with 

P 2 = Ci 2 + C 2 2 , ^tan- 1 ^), (79) 

and 

Ci = pcos#, C 2 = psin6>. (80) 



21 



The eigenfunction Q u j has the form 



(81) 



where T u (p) is given by flT2"|) and 



(82) 



The wave function ty u i e Q,5 2 has the form 



e tn T u>l (p)e 




(83) 



and T u is given by (j72|). 

5.2 Rotation of the Operators 

To obtain the wave function T u> i(p) using the ladder operators a± and c±, we trans- 
form the generators of the os(2) subalgebra into a form that can be written as a 
product of a time-dependent function and an operator depending only on p and 9. 
Let o be some operator which acts on the manifold of solutions Qs 2 - Then, we have 



oe 



e-^^(Ci,C 2 )T(/i) 



e iK e -r,L z e „L z e -iK jn e - V L z ^ , £ )T(//) 



e l7 V^C^(Ci,C 2 )T(//), 



(84) 



where the new operator O is defined to be 



O = e vL *e- in oe in e- vL % 



(85) 
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with 1Z given by ([72]) or ([82]). The operator O acts on the product space of functions 
denoted by {T Ui i(p)e zW T u (p)}. Note that the generator C z of the o(2) algebra is 
invariant under this transformation. 

Let the operator o be a generator of the Heisenberg-Weyl algebra u> 2 . For a±, 
given by (|55|) , we get 

e^e^a-e^' = £(f)~5e w [0 p + p + i<9 e ] 

= A_, (86) 



and 



= A+. (87) 
For the remaining operators c± of the subalgebra, we obtain 

= C, (88) 



and 



e vL. e -<K c+e «K e -riL. = l<t)\ e ld [d p + p-^d e 



= C + . (89) 
In later applications, it is also useful to introduce the operators 

a = l e »[d p + p+id e ] , (J = \e- l9 [-d p + p+ l -d e l 

c = \e- i6 [d p + p- ^de\ , ct = le*[d p + p-ide], (90) 
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where 

A_ = (f)-5a, A + = (f)iat, 

C_ = (|)-5 C) C7 + = (f)5 c t. (91) 

The operators a, a', c, and c* are the analogues of the corresponding operators in 0. 
See the Appendix. 

Finally, the generator m 3 transforms as 

= M 3 . (92) 

Note that all the spatial dependence has been removed in M 3 . It is now a purely 
/^-dependent operator. 

5.3 Reduction to a Radial Equation 

To complete the solution of Eqs. (|T0|), we shall exploit the representation theory of 
the Lie algebra Q = os(2) o o(2). It is more convenient to work in the new basis 

{d 3 ,f 3 ,a±,c±,I} (93) 

where the operators d 3 and f 3 are defined to be 

d 3 = \(m 3 + C z ), f 3 = \(m 3 -L z ). (94) 
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The operators in this algebra satisfy the nonzero commutation relations: 

[m 3 , a±] = ±a±, [a_, a+] = I, 
[m 3 , c±] = ±c±, [c_, c+] = I, 
[£ 2 ,a ± ] = =pa±, [£ z ,e±] = ±c±. (95) 



The detailed representation theory for this Lie algebra is presented elsewhere [pS . 



The specific representation that is of interest to us is denoted | 1 <g> | 1. The 

~2 ~2 

representation space is spanned by a set of vectors {\n, m), n,m G Z + }, where Z + is 
the set of nonnegative integers. The base for the Lie algebra Q acts on this space in 
the following way: 

f 3 \n,m) = (\ + n)\n,m), d 3 \n,m) — (| + m)|n,m), 



a + |n, m) = y/n+ l\n + 1, m), a_|n, m) = y^ln — 1, m), 



c + |n, m) = y/m + l|n, m + 1), c_ |n, m) = Y/mln, m — 1), 

I\n, m) = \n, m). (96) 

The representation space is a normed space and we have 

(n 1 , m'\n, m) = 5 n > n 5 m > m . (97) 

We identify the solution spaces Qs 2 and Qs 2 with the unrotated and rotated repre- 
sentation spaces of the irreducible representation j i <g> j i. 

It is advantageous to transform c/ 3 and f 3 into operators acting on the solution 
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space Qs 2 of the oscillator. We have 

V 3 = e* L *d 3 e-« L * = l{M 3 + /: z ), (98) 
T 3 = e^f 3 e~^ = \{M 3 -C Z ). (99) 

According to fl§6]), the spectra of these operators are 

S(d 3 ) = S(V 3 ) = {m + § : m E Z+} (100) 
S , (/3)=5 r (J 7 3 ) = {n + i:neZ+}. (101) 

Let 0n jm G Qs 2 De the coordinate representation of the vector \n, m) in the represen- 
tation space | i ® T i- Then, we have 

3®n,m (^ "I - 2)®")™-' 

^3©n, m = (m+i)9 n , m . (102) 

We have already seen that the polar coordinate system (p, 9, p) is a natural coor- 
dinate system for this problem. Therefore, we express T> 3 and T 3 in these coordinates: 



^3 = I - \WzP % + l<f3 + ih + 9 fl ], 

Z>3 = fb 3 d M -i<p> 3 p 2 + |^3 + ^Q-^], (103) 



where £^ = — i^. Substituting |103| ) into ( |102j ) and solving the resulting first-order 



partial differential equations, we obtain the result 

6 n , m = e^T n , m (p)e^( m - n )r n , m (/i), (104) 
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where 1Z is given by (|82|) and 

T n , m = ^kl)^ n+m+1) e- iK °. (105) 

The radial functions T ra m (p) are unknown at this point. If we compare the functions 
( |81| ) with ( |104| ) and ( |72|) with ( |105| ), then we see that they are equivalent solutions if 



we make the identification 

u = n + m, I = m — n. (106) 
Note that these solutions are normed. The norm is 

/+oo r+co _ r+oo r+oo _ 

/ dxdy^ niTn (x,y)^f ntm (x,y) = / dCA^o.oCCi, C2)^o,o(Ci, C2), 

-OO J —CO J — OO J ~ CO 

rco rlix _ 

= pdpdOT nm (p)T nm (p). (107) 

Jo Jo 

5.4 Solution of the Radial Equation 

We can use the algebraic structure of Q in the basis (|93| ) and the structure of the irre- 
ducible representation f 1 <g> f 1 to obtain an explicit form for the radial functions. 

_ 2 ~2 

Since the representation f 1 <8> f 1 is bounded below with respect to the spectra 

~2 ~2 

of both /3 and g?3, there exists an extremal state 

*o,o(p,M) = e-^Q Qfi (p,9,p) 

= e-^e^Vo,o(p^)T ,o(/i) 

= e-^^ToMToM. (108) 



Since 



a_^ ,o = c_^o,o = ^ aT , = cT 0j0 = 0, (109) 
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we obtain 

To,o = ^e-f\ (110) 
after normalization. Therefore, the extremal eigenstates are 



u o,o - e 



*o,o = ^e^e-^e-" 2 / 2 . (Ill) 



This means 0o,o = ^o.o- 

According to Eqs. fl96|), we have 



= ^fer<e-^ ,o. (112) 
After some calculation [I9|, the complete wave function can be written as 

V n ,m = e^e-^T nim (p)e 49 ( m - m )T nim ( yU ), (113) 
in which the radial wave function T n rrt (p) is 

T ^( P ) = ^7fe4' m - n|) (p 2 )p |m - n| e- 2/2 , (H4) 

where 

k = |[n + m — \m — n\] (115) 
and the /^-dependent function has the form 

T n M = ^V\i)^™+W e - iK °. (116) 



The last equation may be compared to T u (n) in (|72[). Note that for each n there is 
an infinite degeneracy in m. 
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Recall that m 3 and C z are diagonal in the basis {^ n ,m}, since they are linear 
combinations of d 3 and / 3 . From Q9"3j), we see that 

m 3 = d 3 + f 3 , C z = d 3 -f 3 . (117) 

Therefore, the eigenvalues of m 3 and C z are 

^3*n,m = (W + Tl + l)* n ,m (H8) 

and 

£ 2 ^„, m = (m - 7i)$ n , m , (119) 
where n, m G Z + . This implies that the spectrum of 7713 is 

% 3 ) = {m + n + l:n,m6Z + }, (120) 

where each eigenvalue is (n + m + l)-fold degenerate. The spectrum of L z is 

= {m - n : n, m E Z + } = {0, ±1, ±2, ■ • •}, (121) 

which is as expected since the operator C z generates the group 0(2). Each state can 
also be characterized by the eigenvalues of m 3 and C z . Since each state with eigenvalue 
n + m + 1 is (n + m + l)-fold degenerate, the degenerate states are classified or labeled 
by the eigenvalues of C z . They are 

±(m + n), ±(m + n-2), ±(m + n-4), ± 1 or 0, (122) 

according to whether n + m + 1 is even or odd, respectively. 
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6 Solutions to the Pauli Equation 

At this stage, we can return to the original problem discussed in Sec. 1. The solutions 
of the Schrodinger equation found in Sec. 5 and the symmetry structure developed 
in parallel make possible the solution of the Pauli equation and the identification of 
the super symmetry. In this section, we present the solutions of the Pauli equation in 
a useful form. 

The two-dimensional Pauli equation is reduced in Sec. 1 to Eq. (|10|). 

S 2 ±^±=0 , (123) 

where S 2 ± are given in Eq. flTTl) . The symmetry analysis above for each of S 2 + and 
S2- yields the same time-dependent coefficients £, <pi, (p 2 , and ip^. Multiplying S 2 ± 
by — 1<^3 and using [a_, a + ] = I in Eq. ([64] ) and its h 0+ — > h _ partner gives 

-5^3^ = 2a±a T + £ z -m 3 ±±l, 

= 2a±a T -2/3±±l, (124) 

where m 3 ± and f 3 ± are given by ( p4|) and ( |94"D and their partners. Equation (|l0|) can 
therefore be rewritten as 

(2a±a T - 2/ 3± ± 1)*± = 0. (125) 
From Sec. 5, the solutions to this equation can be taken as 

= e^e-" L ^n ± , m± (Ci,C2)T K±;n± , m± (/i), (126) 
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where ^ n+)m+ (x, y, r) is given by (|113| ) and ^ n _ >m _ (x, y, r) is its partner. The function 
ip n ,m is a function of either of the variables of separation ( 2 ) or (p, 6). The result 
is a double eigenvalue problem: 



a±a T ^ n±tTn± = I ™ + 1 *„ ± , m± . (127) 



n_+l 



Using ( |126| ) for the solution ^ n+tm+ and a similar expression for its partner, and 



noting (^) and (|87D , we obtain the result 



V " • 1 / 

Finally, since the operators A± contain no time derivatives and because of Eqs. 
(PIP, the two-component Pauli equation can be written as 



oat i I ^„_, m _ / I n-+l \ if 



n_ ,m_ 



(129) 

where the operators o)a and aa^ are given in the coordinate representation by 
« f a = lh(5ciCi + ^0i)-2i(CiSci-Cififii) + (Ci + ^)-2], 

= i[-(%c 1 + ^c 2 )-^(C 2 %-Ci%) + (Ci 2 + Cl) + 2]. (130) 

Equation ( |1 29|) is the generalization of the factorization obtained in Ref. . 
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7 Supersymmetry 



In this section, we identify a supersymmetry associated with a nonrelativistic charged 
spin-| particle moving in a time-varying uniform magnetic induction B(r) = B(r)z. 
This supersymmetry generalizes that discussed in Refs. @, ||, [7]]. 

The relevant time-dependent Pauli equation for this system is Eq. ([!]), with 
4>(r,t) = and the gauge choice in Eq. (Q). For more generality, we do not take 
p z as zero but instead separate variables with respect to z. The momentum in the 
z-direction is then represented in the expressions below by its eigenvalue k±. 

As described in the previous sections, the solution space Qs 2 + spanned by the func- 
tions {^ n+t m + : n + ,m + G Z + } forms a basis for a representation space for the irre- 
ducible representation f 1 <E> t I °f the Lie algebra of operators {d 3+ , f 3+ , a±, c±, I}. 
A partner space Qs 2 - spanned by the partner functions {^ r n _, m _ : n-, m_ G Z + } also 
exists and forms a basis for the same irreducible representation of the isomorphic Lie 
algebra of operators {c? 3 _ , / 3 _, a±, c±, I}. The spectra are 

S(fs±) = {n ± + \:n ± e Z+}, S(d 3± ) = {m± + | : m± G Z+}. (131) 

From Eq. ( |129| ) we see that the spectrum of the operator a) a is 

S(a j a) = {n + : n + G Z+} = {0, 1, 2, • • •}, (132) 

while the spectrum of aa) is 

S(aa ] ) = {n_ + 1 : n_ G Z+} = {1, 2, 3 • • •}. (133) 

In effect, the two operators a) a and aa^ have identical spectra, except that the latter 
is missing the ground state. This suggests the existence of a supersymmetry. 
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Following Ref . J7| , we introduce a unified notation that permits the simultaneous 
handling of the two spaces. We define a parameter v that takes the value for the 
'bosonic' space Qs 2+ an d 1 for the 'fermionic' space Qs 2 ~- ft distinguishes the upper 
and lower components of the two-component Pauli equation. States in the two spaces 
can then be denoted by \n,m; is), where n = 0, 1,2, . . . and m + | G Z + . For the 
bosonic space n = n + and m = m + , while for the fermionic space n = n_ and 
m = m_. In two-component notation, we have 

\n, m; v) = ip n , m ((i, C2) J ■ (134) 

The action of the operators a, c and their conjugates becomes 

a\n, m; v) = ^Jn\n — 1, m; u), a^\n, m; v) = + l\n + 1, m; z/), 

c\n,m\v) = ^pm\n,m — \\v) , c' \n, m; v) = + l\n, m + 1; v), (135) 

These expressions are the natural time-varying extensions of the results obtained in 
Ref. 0. 

We can also define raising and lowering operators b and for the index v. By 
definition, these operators satisfy the anticommutation relations 

{b, b ] } = I, {b, b} = {b\ b ] } = 0. (136) 

Their two-component form is 

*-(!!)• (137) 

Their action on the ket |n, m; v) is given by 

b\n, m; v) — 5x„\n, m\ 0), b*\n, m;u) — 5q v \ti, m; 1). (138) 
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We can now introduce the time-dependent operator H, defined by 



H = a ] a + b ] b. 



(139) 



Note that, although H is an integral of the motion, there is implicit time dependence 



in a and a) . When normalized as H = eBH/M, this operator is a time-dependent 



extension of the usual hamiltonian for the time-independent Landau problem. It 
directly generalizes the operator denoted by H in Ref. 0. Using Eq. ( |139|) and the 
new ket notation, the factorized Schrodinger equation ( |129| ) becomes 



This expression shows that the states \n + ,m + ;0) and |n_ = n + — l,m_; 1) are de- 
generate, except for the unique ground state |0,0;0), which has zero eigenvalue for 
H (as required for unbroken supersymmetry) . 

The framework for the supersymmetry is now almost complete. It remains merely 
to introduce appropriate supersymmetry generators mapping degenerate states into 
one another. These operators are defined by 



They satisfy the graded commutation relations of the supersymmetric quantum- 
mechanical algebra sqm(2): 




(140) 



Q = ab\ Q f = a ] b. 



(141) 



{Q,Q} = {Qt ; gt } = , {Q,Q^} = H, 



[H,Q} = [H,Q^] = 0. 



(142) 
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Explicitly, the action of the supersymmetry generators is 

Q\n, m; v) = ^/n5 0u \n - 1, to; 1), 
Q^\n, m\v) = \fn + 15\ u \n + 1, to; 0). (143) 

The operator Q maps bosonic states into fermionic ones, while its conjugate does the 



reverse. 
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8 Supercoherent States 

At this stage, we are in a position to construct the supercoherent states of the time- 
varying Landau system. To do so requires a supersymmetric generalization of the 
standard displacement-operator method. A natural approach to this was introduced 
in Ref. |7j], using the supermanifold formalism developed in Ref. [2^] and the tech- 



niques for BCH relations of Ref. | 23fl . For a description of the general procedure and 



examples of its application, see the above references and Refs. p4| , |25 . 

In the present case, the relevant superalgebra Q s is the one obtained by extending 
the Lie algebra Q as follows: 

Q s = {a f a, c f c, a, a\ c, c\ 1; 6 f 6, b, b ] }. (144) 

A fixed state is required by the construction. We choose it as the ground state |0, 0; 0). 
The subalgebra consisting of the operators {a^a, c^c, 1; tfb} leaves this state fixed, i.e., 
the ground state is an eigenvector of these operators. According to the procedure of 
Ref. J7|, the supercoherent states are to be defined via the action of the operators in 
the quotient algebra {a, a\ c, c^; b, tf} on the fixed state. 

It is convenient to define a super- Hermit ian basis for the quotient algebra: 

Xi = a + aJ, X 2 = i(a — a'), X 3 = c + c* , A 4 = i(c — c'), 

X b = i{b + tf), X 6 = tf-b. (145) 

Constructing a unitary representation of the supergroup for an element in the quotient 
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algebra |7|, |25|, we obtain 

T(g) = exp(iA l X 1 +iA 2 X 2 + iC 3 X 3 + iC 4 X A + i9 1 X 5 +i9 2 X 6 ), 

= exp(-Aa + Aa j -Cc + Cc j + 6tf + 6b), (146) 

where A = A 2 + iA x , C = C 2 + id, and 9 = -61 + i6 2 . This means that A, C G °B L 
and 9 G 1 Bl are complex Grassmann- valued variables. 

Using a suitable Baker-Campbell-Hausdorff relation for the supergroup element 



and Lemma 1 of Ref. |[23|| , we can construct the analogues of the supercoherent states 
in Ref. |7|]. The states are parametrized by three Grassmann- valued parameters A, 
C, and 9, and are given as 

\Z) = exp(±^)exp(-±|,4| 2 )exp(-±|Cf) 



Recall that the operators a, cr, c, and have an implicit time-dependence in them 
through their dependence upon the variables Ci = X / ( P^ 2 ( T ) an d C2 = y/ l pl^ 2 ( T )- The 
supercoherent states defined above are the natural time- dependent generalizations of 
the supercoherent states derived in J7|. 

In Ref. |26|| , several theorems about time- dependent integrals of the motion are 
given. The first states that any function of integrals of the motion is itself an integral 
of the motion. The second states that eigenvalues of time-dependent integrals of the 
motion do not depend on time. The third states that applying an integral of the mo- 
tion on a solution to a wave equation (either a Schrodinger or a Pauli time-dependent 
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equation) yields a function that itself is a solution to the same wave equation. Since 
a, b, c, and their conjugates are integrals of the motion, so too is the operator T(g). 
Furthermore, since |0, 0, 0) is a solution to the Pauli equation, \Z) is also a solution to 
the same equation. All the formulae in Ref. may be rewritten identically for the 
properties of the present supercoherent states, with obvious replacements. Note that 
some aspects of integrals of motion in simple problems with supersymmetric quantum 



mechanics have been discussed in Ref. [27 



If one introduces a time- dependent electric field in addition to the present mag- 
netic field, the basic algebraic structure derived here can again be applied. Provided 
this leads to a supersymmetric formulation, this provides a means of obtaining the 
explicit form of supercoherent states for this more general case. It is plausible that 
supercoherent states for a relativistic electron moving in an electromagnetic field of 
arbitrary configuration could be found using these methods combined with coherent 



states in the proper-time formalism, as discussed in Ref. [2g . 
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9 Appendix: Constant Magnetic Induction 

In this appendix, we demonstrate explicitly the connection between our results and 
the special case of constant induction treated in Ref. (Tj. 

When the magnetic induction B is constant and uniform, w = eB, the differential 
equations ( p6|) and (|27|) take the form 



e Q + \w 2 e a = 0, a = 1, 2, (148) 

and have real solutions 



Xi = V5 X2 = v / |^( 2 f) (149) 



satisfying the wronskian condition W(xi) X2) = 1- The complex solutions ( fnj ) become 
e = 75(Xi + %) = ^e-/ 2 , £ = ^~ lWTl \ (150) 

with wronskian = — z. 

The complex solutions to the equation ([25]) are given by 



¥>i = £e*"\ ^ 2 = ^e— , ^3 = ^, 17 =|«w. (151) 



Using the equations derived in Sec. 5.2, we can write down the generators of the 
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Schrodinger algebra 

m~ = ±e™ T [-d T -^L z -f(xd x + yd y )-iw 2 (x 2 + y 2 )-f-ih Q ], 
m + = l e -™ T [9 T + fL 2 -f ( X d x + yd y ) + iw 2 (x 2 + y 2 )-f + ih ], 
™3 = %[d T + iL x + iho], 



C z = iL z = i(yd x - xd, ; 



yji 



a - = b^S^x + i9 v + + 



a + = ^[-( 9 x ~ id y ) + ^(x - iy)}, 



c + = -^[-(d x + id y ) + ^(x + iy)), 
1 = 1. (152) 
The generators of the subalgebra Q are Q = {d 3 , fo, a±, c±, I}. We have 

/a = ~(i8r - h ), 

d 3 = ±(id T + w£ z - h ), (153) 

along with the operators (|152|) . 

From the results of Section 5, the solutions ( |1 13| ) are 

V n , m = e-^T n>m (p)e id( - m ^T n>m (fx), (154) 

where T n>m (p) is given by (|114j) and 

T = (— )5 e - i ( n+m+l ) w t 1 / 2 e - ih ot j - (155) 
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where n, m £ Z + . The polar coordinate system is related to the cartesian system 
through the transformation 

P 2 = Ci+Cl, = tan- 1 C 2 /C2, (156) 

where 

Ci =x(f)l, Ci = x(f)i, At = r , (157) 

following Eqs. (|TT]) . 

By direct calculation, we see that the action of f% and d 3 on \I/ nim is 

d 3 ^ n>m = {m + \)-$ n>m . (158) 

Also, the ladder operators a± and c± raise and lower the quantum numbers n and m, 
respectively. Thus 



In the special case of constant field, the r or p variable may be separated and 
we can work with the cartesian-like coordinates (£i, (2) or with the polar coordinates 
(p, 9), as defined in (|157|) and (|156|), respectively. In the former case, use the wave- 



function ijJn,m{Cii C2)) an d in the latter we use ip n ,m(p, 0) = T njm (p)e* e ( m n \ On these 
solution spaces, the operators / 3 and d 3 are replaced by N a = a^a and N c = c^c, 
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respectively. The ladder operators a and c are defined by 

a = |[(9 Cl +z%) + (Ci + <2)], 

c = |[(%-^ C2 ) + (Ci-< 2 )], 

= vfe[(^-^) + f(^-^)], (160) 



where we have made use of ( |157| ). In polar coordinates (|156|) , the operators and their 
conjugates are given in Eqs. (|90D . 



The operators in ( |160| ) are exactly the a and c operators obtained in Ref. 0. 
Their time-dependent extensions are a_ and c_, respectively. 

The number operators obey the nonzero commutation relations 

[N a , a] = -a, [N a , at] = +a\ [a, a)} = I, 
[N c ,c] = -c, [JV c ,ct] = + C t, [c,ct]=J. (161) 

The action of these operators on the manifold of states {ipn,m} is 

N a 1p n , m = nijj n , m} aipn :m = V^n-l.m, otyn.m = V n + fyn+l,m, 

Substituting ( |160| ) for a in N a , we obtain 

Na = i[-(9 ClCl + (% C2 +2^) + (Ci 2 + C 2 2 )-2] 

= ±[-(d xx + d yy + wC z ) + ^(x 2 + y 2 )-w], (163) 
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which is proportional to the Hamiltonian |7]]. 

This completes the explicit demonstration that our formalism contains as a special 
case the description of the motion of a nonrelativistic spin-| charged particle in a 
constant and uniform magnetic induction. 
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